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ABSTRACT 

We numerically explore planet formation around a Centauri A by focusing on the crucial 
planetesimals- to-embryos phase. Our a pproach is significantly improved with respect to the 
earlier work of Mar zari & Scholll d2000l) . since our deterministic N-body code computing the 
relative velocities between test planetesimals handles bodies with different size. Due to this 
step up, we can derive the accretion vs. fragmentation trend of a planetesimal population 
having any given size distribution. This is a critical aspect of planet formation in binaries 
since the pericenter alignment of planetesimal orbits due to the gravitation al perturbations of 
the co mpanion star and to gas friction strongly depends on size. Contrary to Marzar i~& Scholll 
(2000), we find that, for the nominal case of a Minimum Mass Solar Nebula gas disc, the 
region beyond ~ 0.5 AU from the primary is strongly hostile to planetesimal accretion. In this 
area, impact velocities between different-size bodies are increased, by the differential orbital 
phasing, to values too high to allow mutual accretion. For any realistic size distribution for 
the planetesimal population, this accretion-inhibiting effect is the dominant collision outcome 
and the accretion process is halted. Results are relatively robust with respect to the profile and 
density of the gas disc. Except for an unrealistic almost gas-free case, the inner "accretion- 
safe" area never extends beyond 0.75AU. We conclude that planet formation is very difficult 
in the terrestrial region around a Centauri A, unless it started from fast-formed very large 
(>30km) planetesimals. Notwithstanding these unlikely initial conditions, the only possible 
explanation for the presence of planets around 1 AU from the star would be the hypothetical 
outward migration of planets formed closer to the star or a different orbital configuration in 
the binary's early history. Our conclusions differ from those of several studies focusing on the 
later embryos-to-planets stage, confirming that the planetesimals-to-embryos phase is more 
affected by binary perturbations. 

Key words: planetary systems: formation - stars: individual: a Centauri- planets and satel- 
lites: formation. 



1 INTRODUCTION 

A majority of stars are memb ers of binary or multiple systems (e.g. 
Duquennov and Mavor] [l99 ll) . making the presence and frequency 
of planets in such systems, in particular terrestrial planets, a very 
relevant issue in our ongoing search for other worlds. This issue 
has been made all the more timely by the discovery of more than 
40 extrasolar planets in multiple systems, amounting to ~ 20% of 
all de tected exoplanets (Rag havan et al .12006k Desidera XBarbieril 
120071) . Most of these planets are located in very wide binaries 
(> 100AU) where the companion star's influence at the planet's lo- 
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cation is probably weak (except for possible high incli nation effects 
betwe en the planet's and the binary's orbits, see e.g. Tak eda et al.l 
120081) . However, a handfull of planets, like Gliese 86, HD41004 or 
y-Cephei, inhabit binaries with separations ~20 AU, for which the 
binarity of the system probably cannot be ignored. 

The long term stability of planet ary orbits in binaries has 
been explored in several studies (e.g. iHolman & Wiegertl [l 999 ; 
David et al. I l2003l ; iMudrvk & Wi] l2006). aimed at estimating a„ it , 
the distance from the primary beyond which orbits become unsta- 
ble. This critical distance depends on several parameters, mainly 
the binary's semi-major axis a b and eccentricity e b as well as the 
mass ratio fi = m 2 /(' )I i + m i)- However, one can use as a rule of 
thumb that, except for extremely eccentric systems, a all is in the 
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0.1-0.4flj, range (see for e xample the empirical formulae derived by 
IHolman & Wiegert| [T999). As a consequence, the terrestrial planet 
region around 1 AU is almost always "safe" for binaries of separa- 
tions > 10AU. 

The way the planet formation process is affected by stellar 
binarity is a more difficult issue. T his process is in deed a complex 
succession of different stages n(e.g |Lissauelll993h . each of which 
could react in a different way to the companion's star perturbations. 

The last stage of planet formation, leading from plan- 
etary embry os to planets , has been investigat e d in several 
papers (e.g IBarbieri et all 120021 : lOuintana et ail |2002| . 120071 ; 
Haghighipour & Raymond 2007). These studies have revealed that 
the (a b ,e b ,(i) parameter space for which this stage can proceed 
unimpeded is slightly, but not too much more limited than that 
fo r orbital stabili t y. As an illustration, the numerical exploration 
of lOuintanaetal J ( 120071) showed that planetary embryos can mutu- 
ally accrete at 1 AU from the primary if the companion's periastron 
q2is> 5 AU. 

In contrast, the stage preceding this final step, i.e. the mu- 
tual accretion of kilometre- sized planetesimals leading to the em- 
bryos themselves, is much more sensitive to binari t y effects. A se- 
ries of numerical studies llMarzari & Scholl 2000; Thebau lt et al.l 
l2004l2006l : |Paardekooper et alj2008t) have indeed shown that per- 
turbations by the companion star can stir up encounter velocities 
<Av) to values too high for small planetesimals to have accret- 
ing impacts. Such accretion-inhibiting configurations are obtained 
for a wider range of binary para meters than for the lat er embryo- 
to-planet phase. As an example, iThebault et alj J2006t) (hereafter 
TMS06) have shown that planetesimal accretion in the 1 AU region 
can be perturbed in binaries with separations as large as 40 AU. 
These results are not surprising in themselves, and this for two rea- 
sons. 1) The growth mode in the planetesimal-accumulation phase 
is the so-called runaway growth, whose efficiency is very sensitive 
to increases of the encounter velocities (which reduce the gravi- 
tational focusing factor, e.g. TMS06). 2) It takes a much smaller 
perturbation to stop accretion of a 1km body than of a lunar-sized 
embryo: (Av) >10m.s~' is typically enough to achieve the former, 
while large embryos will hardly feel any difference in their accre- 
tion rate for such small velocity disturbances. A crucial mechanism 
driving the (Av) evolution of the planetesimals is friction due to the 
primordial gas remaining in the protoplanetary disc. In a circumpri- 
mary disc, its main effect is to induce a strong phasing of neighbor- 
oughing planetesimal orbits, which cancels out the large orbital pa- 
rameters oscillations around th e forced eccentricity ef induced by 
the companion's perturbations dMarzari & S choll 2000). A crucial 
point is that this effect is size dependent, so that while (Av) between 
equal-sized objects are reduced, impact velocities between objects 
of different sizes are increased. TMS06 have shown that the dif- 
ferential phasing effect appears already for small size differences, 
suggesting that this effect may dominate in "real" planetesimal sys- 
tems. However, the exact balance between this accretion-inhibiting 
effect and the accretion-friendly effect on equal-sized objects de- 
pends on the size distribution in the initial planetesimal popula- 
tion, a parameter which could not be thoroughly investigated in the 
global investigation of TMS06. 



1.1 planet formation in a Centauri A 

We follow here the approach of TMS06 to investigate in detail one 
specific system: a Centauri A. Focusing on one single system al- 
lows to take the TMS06 model a step further and explore realistic 
size distributions of the planetesimal population as well as crucial 
parameters such as gas disc density and radial profile. 

At a distance of 1.33 pc, the a Cen AB binary sys- 
tem, with a possible third d istant companion M-star Proxima 
dWertheimer & Laughli d2006T) . is our closest neighbour in space, a 
Centauri A and B are G2V and K1V stars, of masses M A = 1.1M Q 
and M B = 0.93M o respectively. The binar y's semi-major axis i s 
a b = 23.4 AU and i t s eccen tricity e b = 0.52 jPourbais etal]|2002l) . 
IHolman & Wiegertl d 1997b have shown that the region within ~ 
3 AU from the primary can harbour stable planetary orbits. How- 
ever, the search for planetary compa nion has so far b een unsuc- 
cessful. The radial velocity analysis o fiEndl et al.l d200lh has shown 
that the upper limit for a planet (on a circular orbit) is 2.5 M ]upiler 
around a Cen A, and this at any radius from the star. This leaves an 
open possibility for the presence of terrestrial planets in the crucial 
~ 1 AU region. 

Planet formation ar ound a Cen A has so far b een specifically 
investigated in 4 studies. IBarbieri et al.l d2002h and lOuintana et all 
(2002, 2007) have focused on the late embryos-to-planets stage and 
have found that it can proceed unimpeded in a wi de inner region (< 
2.5 AU ). For the planetesimal-accretion phase, iMarzari & Scholll 
(2000) have found that orbital phasing induced by gas drag allowed 
km-sized bodies to accrete in the r < 2 AU region. However, this 
study only considered impacts between equal-sized objects, thus 
implicitly overlooking the accretion-inhibiting differential phasing 
effect. We reinvestigate here this crucial phase, exploring in detail 
gas drag effects on realistic size distributions for the planetesimal 
population. 



2 MODEL 
2.1 algorithm 

We use the same deterministic N-body code as in TMS06, which 
follows the evolution, in 3D, of a swarm of test particles submitted 
to the gravitational influence of both stars and to gaseous friction. 
This algorithm tracks all mutual encounters within the swarm and 
allows to derive precise statistics of the evolution of the (Av) distri- 
bution. Gas drag is modelled by assuming a non evolving axisym- 
metric gas disc, where the only free parameter is p g (r), the radial 
profile of the gas density distribution. The drag force then reads 



F 



3p g C d 



-vv, 



(1) 



where v is the velocity of the particle with respect to the gas 
and v the velocity modulus. p pl and s are the planetesimal den- 
sity and physical radius, respectively. Cd is a dimensionless co- 
efficient = 0.4 for spherical bodies. Note that the planetesimal's 
physical size .? considered for the drag computation cannot be taken 
for the close encounter search routine, because the collision statis- 
tics would be much too poor with the limited number of particles 
(typically 10 4 ) inherent to N-body codes. We thus have to rely 
on the usual prescription of assigning an inflat ed radius s num to 
the particles for the encounter search routine < Thebault & Brahij 
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20071) . This proce- 



1 we consider here the widely accepted "core-accretion" scenario 



dure is valid as long as it does not introduce a bias, in the form of 
an artificial shear component, in the (Av) statistics. We take here 
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10 5 AU, a value which gives a precision of the order of 



1 m.s 1 in the encounter velocity estimate^. 



2.2 gas disc prescription 

The axisymmetric prescription is a simplification of the real be- 
haviour of the gas disc, which should also react to the companion 
star's perturbations and departs from perfectly circular streamlines 
dArtvmowicz and Lubow|[l994l) . To realistically model the evolu- 
tion of both the planetesimal and gaseous components, using a cou- 
pled N-body and hydro approach, is an ardu ous task. First attempts 
at addressing this issue have been done by Ciecielag et alj l2007h 
for a simplified circular-binary case and b y Paardekoop er et all 
(2008) for eccentric binaries. IPaardekooper et al d2008l) have 
shown that the gas disc becomes eccentric under the companion's 
perturbations but that the amplitude of this eccentricity e g depends 
on the choice for the numerical wave damping procedure (the flux- 
limiter prescription). There are two possible modes: a "quiet" low- 
e g case for first order flux limiters and an "excited" high-e s case 
for second order ones. In the "quiet" regime, the (Av) statistics 
for planetesimals is relatively close to the simplified axisymmet- 
ric case. In fact for a binary of 20 AU separation, results between 
the two cases were indistinguishable. In the "excited" regime, en- 
counter velocities are on average increased by a factor of ~ 2, 
with respect to the circular case, for a 20 AU separation binary. It 
is important to point out that, regardless of the details of the wave- 
damping modelling procedure, the (Av) are either equal to or higher 
than their v alues in the circular disc c ase (for a more detailed dis- 
cussion, see Paardekooper e t alfcooih . 

Because of these uncertainties regarding the "real" behaviour 
of the gas disc, we chose to keep here the simplified axisymmet- 
ric gas disc approach. This choice is also motivated by the fact 
that these full planetesimal+gas simulations are still very CPU ex- 
pensive and are thus not suited for detailed parameter-space explo- 
rations as the one which is t he purpose of this study. H owever, in 
the light of the conclusions o f lPaardekooperet al. (2008), we think 
that our results should not be too different (within a factor 2 in 
terms of encounter velocity values) from those expected with a full 
N-body + Hydro modelling. In any case, a crucial point is that our 
results should be regarded as giving a lower boundary for the (Av) 
estimates. 



Table 1. Nominal case setup. The fields marked by a V are explored as free 
parameters in the simulations. See text for details on each parameter. 



2.3 set up 

We consider for each run a number N = 10 4 of test particles, each 
having the same numerical size (inflated radius) s num = 10~ 5 AU 
and physical sizes (for the gas drag computation) taken at random 
between two boundaries s mj „ and s max . For our nominal case, we 
consider planetesimals in the s mj „ = 1 km to s„ mx = 10 km range. 
We also consider an alternative "small planetesimals" case with 
s min = 0. 1 km and s max = 1 km and a "large planetesimals" case 
with s mjn = 5 km and s max = 50 km (see Sec |3.4| >. The test par- 
ticles are initially distributed in an approximate terrestrial planet 
region, between r mj „ = 0.3 AU and r„„-„ = 1.5 AU. They are ini- 
tially placed on almost circular orbits ((e)o < 10~ 5 ), so that initial 
encounter velocities are low enough to be in the accretion-friendly 
range. This choice might not correspond to a physically realistic 

2 to a first order, the error in the (Av) estimate is of the order of v/c, p (r) X 
s„ um /r, where r is the radial distance to the star and \>Kep the orbital velocity 
at distance r 
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0.3-1.5 AU 




Number of test particles 


10 4 


V 


Physical sizes s mm ,s max 


l-10km 




Inflated radius (collision search) 


irr 5 AU 


V 


Gas disc radial profile 


Pg =p g0 (r/lAUr 2 - 75 


V 


PgO 


1.4 x KT'g.cnT 3 




ffinal 


10 4 yrs 



case, but TMS06 have shown that the initial orbital conditions are 
very quickly relaxed and that the system tends towards the same 
steady-state regardless of the initial (e) and (/). 

For the gas density, we consider as a nominal case the canon- 
ical Min imum Mass So lar Nebulae (MMSN) radial profile de- 
rived by iHavashil dl98ll) . with p g = p g0 (r/lAU)" 2 - 75 and p g0 = 
1.4x 10~ 9 g.cirT 3 , but explore different radial profiles and p g o values 
(see Sec |3.3| l. 

The final timescale for each run is taken to be fg„ = 10 4 yrs, 
i.e. the typical timescale for the runawa y accretion of kilometre- 
sized planetesimals (e.g. lLissauedfl993l) . However, this choice is 
not too stringent, since in most cases a steady-state (Av) distribution 
is reached in less than a few 10 3 yrs. 

All characteristics of the nominal case set-up are summarised 
in TabQ 



2.4 (Av) analysis 

The distribution of encounter velocities is the main output we are 
interested in. For each mutual impact found by the close-encounter 
search routine, we record Av sl>S2 , si and S2 being the (physical) 
sizes of the impacting bodies. The encounter velocity distribu- 
tions are then classified in 2-D matrixes giving Av averaged over 
[si, si + ds] and [s2, si + ds] size intervals in the [j„„„, s max ] range. 

The crucial issue is of course to interpret these velocities in 
terms of accreting or eroding impacts. For an impacting pair of 
sizes and s%, the limit between eroding and accreting impacts can 
be defined by the threshold velocity v*( SljS2 ). This velocity depends 
on several parameters obtained from energy scaling considerations. 
The main one is the threshold energy for catastrophic fragmentation 
<2*B for which vari ous and often conflicting esti mates can be found 
in the literature (see Theba ult & Augereaul2007l for a thorough dis- 
cussion on the subject). We follow the same conservative approach 
as in TMS06 and consider 2 boundary estimates for the threshold 
velocity, v* (j] s , )low and v* (J1 S2)high , derived using a weakly-resistant 
and strongly-resistant Q* prescription respectively (see TMS06 for 
more detail). For velocities in between these 2 values, we do not 
draw any conclusions regarding the accretion vs. erosion balance. 

Note that, despite their diversity, all prescriptions for Q* agree 
that the size ranges considered here, i.e. 0.1 to 50 km, correspond 
to the so-called gravitational regime, where bodies' resistance is 
determined by gravitational binding forces. This means that Q* is 
an increasing function of size. This point is crucial for the results 
of Sec|3~4l 



3 although it is important to stress that cratering (i.e., non fragmenting) 
impacts play a crucial role in defining the accretion/erosion limit, the pa- 
rameters for cratering, mainly the excavation coefficient, follow (to a first 
order) similar energy scaling laws as Q* 
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Figure 1. Mean encounter velocities (Av si , s ,), averaged in the 0.9-1.1 AU 
region, for all impacting pairs of sizes S[ and s 2 in the [1, 10km] range. For 
the sake of clarity, 4 iso-velocity lines are represented. 
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Figure 2. Same results as in Fig[T] but to each (Av S[ , S2 ) value is associ- 
ated a colour scale, depending on how this value translates into accreting or 
erosive impacts for each impacting pair of sizes si and s%: 
-green: no (Av) variation observed, region of unimpeded accretion 
-blue: encounter velocity increase but (Av) < v*i ow , perturbed accretion 
with slow-down of the runaway growth 

-yellow: v*i ow < (Av) < v* n igh. uncertain accretion vs. erosion balance 
-red: (Av) > v*^, erosion 
-black: (Av) > 2 X v*high> erosion 

3 RESULTS 

3.1 (Av) distribution, erosion threshold 

FigUshows the (Av S| , S2 ) distribution, at 1 AU from the primary and 
for bodies in the [1, 10km] range, once a steady-state is reached, 
i.e. after = 3 x 10 3 yrs. The sharp size dependence of encounter 
velocities is here clearly visible. Low (Av) are restricted to a nar- 
row stripe along the s\ = S2 diagonal, corresponding to the domain 
where gas drag aligns orbits towards the same value. However, for 
any small departure from the exact s t = s 2 condition, encounter ve- 
locities dramatically increase. We get (AVji,^) > 50m.s~' in almost 



Figure 3. Relative importance of different types of collision outcomes (cor- 
responding to the different colour categories 0fFig[2}, as a function of dis- 
tance to the primary star. All close encounters between all bodies in the 
[Smin> s max] s i ze range are taken into account, with each close encounter 
weighted assuming that the size distribution for the planetesimals follows a 
Maxwellian distribution centered on 5km (see text for details). 

the whole s\ t S2 region. This illustrates very clearly the strength 
of the differential orbital phasing mechanism. 

FigO shows how these velocities can be interpreted in terms 
of accreting/erosive impacts. The only region where encounter ve- 
locities are low enough for the accretion process to proceed unim- 
peded (in green) is the narrow si = S2 diagonal. This diagonal 
is surrounded, for objects >5km, by an equally narrow band (in 
light blue) where (Av) allow accretion ((Av) s , s , ^ v* (n _ r2) i ow ) but 
are nevertheless increased with respect to their initial unperturbed 
value. This means that the accretion rate is slowed down, com- 
pared to the unperturbed case, since the runaway growth mode's 
efficiency decreases with increasing encounter velocities. Except 
for a small region (in yellow) of uncertain accretion/erosion bal- 
ance where v* (JM2)low < (Av) S|>s , < v* (J1 , 2)high , all other (s u s 2 ) 
encounters result in eroding impacts. In these red and black areas, 
which cover most of the graph, velocities exceed v*( sllJ2 )iiigh> our 
most conservative estimate for the threshold velocity . 

3.2 accreting vs eroding impacts: size distribution 

Of course, Fig(2]cannot in itself give the relative importance of the 
accreting-friendly vs. accretion-hostile areas. To do so, the real size 
distribution of the planetesimals should be known, whereas in the 
runs all sizes were chosen at random between s mj „ and s„ mx . The 
size distribution of a typical initial planetesimal population is a 
complex issue, which is directly linked to the very difficult (and 
yet unresolved) problem of the mechanisms for planetesimal for- 
mation. This problem exceeds by far the scope of the present paper 
(see, however, the discussion in TMS06 and in the last Section of 
the present paper). In line with our conservative approach, we shall 
here assume a fiducial size distribution which is a priori favourable 
to same-size impacts (that is, to accretion), i.e. a Maxwellian f mx (s) 
peaking at a median value s med = 5 km. 0. We then weight each 

4 Some planetesimal-accumulation particle-in-a-bo x models make the as- 
sumption of a single-size initial population (e.g. IWeidenschilli ntTet alj 
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(si,i2) impact by fi = f mx (si) x f mx {Si), the (normalised) probabil- 
ity of such an impact with a f mx {s) size distribution. Depending on 
<Av) S j _ s , , all impacts are then classified into one of the 5 collision- 
outcome categories defined earlier (see FigO, where their contri- 
bution is weighed by fi. Comparing the total E/j for each category 
then gives the relative importance of each of these possible colli- 
sion outcomes. 

As can be seen in Fig[3] at 1 AU less than 2% of all impacts 
result in accretion (a value which only increases to ~ 5% when in- 
cluding the "uncertain-outcome" category), while 95% should lead 
to net erosion of the impacting bodies. The dynamical environ- 
ment is thus globally very hostile to accretion. Fig[3]also shows that 
this accretion-hostile environment extends over most of the system. 
Only the innermost r < 0.5 AU region has a favourable accretion 
vs. erosion balance. This accretion-friendly inner zone is mostly 
due to weaker gravitational perturbations by the companion star. 

Between 1.2 and 1.4AU, a slight increase in the accret- 
ing/eroding impacts ratio is observed, which is mostly due to 
weaker gas drag-induced differential phasing. However, even there 
the fraction of impacts with (Av) < v*i ow never exceeds 8%. More- 
over, beyond 1 .4 AU, accreting impacts quickly decrease to again. 
The reason for high (Av) in these outermost regions is no longer 
gas friction but pure gravitational effects. Indeed, although in the 
absence of gas drag the companion star forces secular orbital os- 
cillations around ey which should in principle not induce high (Av) 
because they preserve the same phasing for all bodies (regardless 
of their size), these oscillations get narrower with time and at some 
point orbital crossing occurs, leading to an abrupt relative velocity 
increase (for a more detailed description of this effect, see section 2 
of TMS06). Strong gas drag can prevent this crossing from happen- 
ing, by cancelling out the oscillations and inducing it's own size- 
dependent phasing, but in the outer regions it is too weak to do so. 
In other words, planetesimals are here caught between a rock and a 
hard place: 

• In the outer regions, pure gravitational perturbations lead to 
orbital crossing and high (Av) 

• In the inner regions, gas drag can prevent orbital crossing but 
induces size-dependent orbital phasing which also leads to high 
(Av) 

We have tested the robustness of our results with respect to 
different choices of the size distribution. Collisional-"equilibrium" 
power laws with dN cc s~ 33 ds and Gaussians centred on s me j of var- 
ious widths have been considered. Although some differences were 
observed with respect to our nominal case, they remain relatively 
marginal and our main conclusions seem to hold regardless of the 
size distribution (see Fig[4^. This is true in particular for the region 
around 1 AU which is always strongly hostile to accretion: the frac- 
tion of accreting impacts never exceeds 10% for all tested size dis- 
tributions. The only case giving an accretion-friendly environment 
is an almost Dirac-like Gaussian of variance (0.25 km) 2 , but this 
size distribution is probably highly unrealistic. Another important 
robust result is that there is always an inner region having low (Av) 
enabling accretion. The outer edge of this inner zone varies slightly 
with the size distribution prescription, but is always located in the 
0.5 -0.6 AU range. 



fl997h , but this is of course a numerically convenient simplification which 
cannot correspond to a "real" population which necessarily has a spread in 
sizes 




Distance to the primary (AU) 

Figure 4. Fraction of accreting impacts (limit between the "blue" and "yel- 
low" zones) for different prescriptions for the size distribution. The noisi- 
ness of the curve corresponding to the extremely narrow-Gaussian is due to 
the fact that in this case only impacts for a very narrow s\ ~ sj region have 
non-negligible weigh in the global statistics. 
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Figure 5. Fraction of accreting impacts (limit between the "blue" and "yel- 
low" zones) for different prescriptions for the gas disc (see text for details). 



3.3 Different gas disc profiles 

For our nominal case, we assumed a standard MMSN gas disc, but 
other gas disc models have been explored (Fig|5j- 

We first consider the same p e = p g o(r/lAU)~ 2 75 profile as in 
the nominal case but with a different reference value for p g0 at 1 AU 
chosen around the MMSN p g o( re fi = 1-4 x 10~ 9 g.ctrT 3 value. 

As it can be seen, for larger values of p 5 o the situation be- 
comes slightly more favourable to accretion, and the outer edge 
of the inner accretion-friendly zone moves to ~ 0.75 AU for the 
highest p g o = 10 x PgO(ref) explored value. This is due to the fact 
that complete differential phasing, meaning that objects of all sizes 
reach steady-state eccentricity and periastron longitude a) values 
before ? nna i, is already achieved for the nominal MMSN case. Thus, 
when increasing the gas density, no "better" phasing is achieved. 
The effect of increased p g is then that all steady-state e and a> val- 
ues globally tend towards their asymptotic value for infinitely small 
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Figure 6. Same as Fig[2] but for the "small planetesimals" run. The colour 
scaling for the different accretion/erosion zones is the same as in Fig[2] 
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Figure 7. Same as Fig|2] but for the "large planetesimals" run. The colour 
scaling for the different accretion/erosion zones is the same as in Fig|2] 



sizes, i.e. e = and oj = 270" (see Eqs. 24-29 of Paardekoop er et al.l 
2008), which tend to reduce encounter velocities. 

The effect of decreasing gas densities is more complex to in- 
terpret. For p g o just below p g o( re f), relatively good differential phas- 
ing is maintained but all steady-state e and oj values globally move 
away from their asymptotic s — > value. This results in a wider 
spread of steady-state orbital parameters and thus higher (Av). Be- 
low a threshold value (~ 0. lxMMSN), however, steady state phas- 
ing is no longer achieved for the biggest bodies, which then still 
have substantial e and u> oscillations after tf, m \. Since phasing is 
still achieved for the smallest objects, the difference between the 
orbits of small and big bodies is still increasing. This is what is 
clearly observed for the p g o = 0. lp g o(ref) run. However, below a sec- 
ond threshold value (~ O.OlxMMSN), steady state phasing breaks 
down for all bodies. The system now moves towards its gas-free be- 
haviour, with large common orbital oscillations for all bodies, and 
high (Av) only in the regions where oscillations are narrow enough 
for orbital crossing to occur, i.e. beyond ~ 1.3AU (see discussion 
in the previous section). 

We finally consider a case with a different radial dependence 
for the gas density, i.e. a flatter profile in p g = p g0 (r/lAU) -1 - 75 
(p g o being the same as in the nominal case). The accretion/erosion 
balance is slightly less favourable than in the nominal MMSN case, 
with the outer edge of the inner "safe" zone at ~ 0.45 AU (dashed 
curved in FigO. 

However, despite all these quantitative differences, the qual- 
itative results remain relatively similar for all test gas disc cases, 
with the exception of the probably unrealistic very-gas-poor runs. 
The generic characteristics are an inner accretion-friendly zone, 
whose outer edge varies from 0.35 to 0.75 AU, followed by a vast 
region where impacts are predominantly erosive for all objects in 
the 1-10 km range. More importantly, the crucial region for hab- 
itable planets around 1 AU is always in an environment hostile to 
accretion. 



3.4 Alternative planetesimal size-ranges 

We have so far assumed a planetesimal population in the 1-10 km 
size range, which corresponds to the most generic choice for an 
"initial" planetesimal population (see e.g. Greenberg et al.|[l978l : 



I Wetherill & Steward ll989UBarge & Peliatl 1 19931 ; llissauerl 1 19931: 
IWeidenschilling et alJI 1997h . However, because the planetesimal- 
formation phase is sti ll po orly understood (see fo r exam ple 
lYoudin & Chiang! J2004I) and lCuzzi & Weidenschillingl j2006h for 
conflicting views on the matter), this few-km size range is far from 
being a value carved in marble. Two other possible initial size 
ranges, with 0.1 < s < 1 km and 5 < s < 50 km, have thus been 
explored in our simulations. Note that if from a purely dynami- 
cal point of view changing the planetesimal sizes is equivalent to 
changing the gas density (because of the l/s dependence of the gas 
drag force), it is not equivalent from the point of view of the accre- 
tion/erosion behaviour. As an example, an impact between a 1km 
and a 5 km body in a gas disc of density p g occurs on average at the 
same velocity (Av)i 2 as an impact between a 10km and 50km body 
in a 10 x p g disc, but the outcome of a (Av)i_2 collision between 
lkm and 5km objects will not be the same as for 10km and 50km 
objects 

For our "small" planetesimals population in the 0.1-1 km 
range, we obtain, not surprisingly, a situation which is more hos- 
tile to accretion (Fig[6j. This is not due to higher relative velocities 
than in the nominal case, since (Av)i^ are of the same order than 
in the 1-10 km run. The explanation is to be found in the fact that 
such small objects are much more affected by (Av) increase, and 
this for 2 reasons: 1) They have lower threshold Q* values, and 2) 
they have lower escape velocities. 

For the "large" planetesimals run with s mi „ = 5 km and 
Simu = 50 km, on the contrary, we find an environment which is 
more favourable to accretion than in the nominal case (Fig[7j>. The 
accretion-friendly zone now covers a large fraction of the (si, si) 
phase space at 1 AU, especially in the > 30 km range. As for the 
small planetesimals case, this result is only marginally due to lower 
relative velocities, since the distribution of (Av)i,2 is here again 
close to that of the nominal case. What makes a difference here 
is that bodies in the 10-50 km range are much more resistant to im- 
pacts than smaller ones. They are also big enough to reaccumulate 
a substantial fraction of the mass they might lose after an impact 
(a process which is taken into account in our estimate of v*( J1(S2 )i 0W 
and v* (J1-S2) hi g h)- 
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4 DISCUSSION AND CONCLUSIONS 

Our numerical exploration has shown that for our nominal case set- 
up the r > 0.5 AU region around a Centauri is hostile to kilometre- 
sized planetesimal accretion. The coupling of the companion star's 
secular perturbations with the size-dependent phasing imposed by 
gas drag leads to high (Av) (si>S2) destructive collisions for all Si + s% 
pairs. For any realistic size distribution for the planetesimal popula- 
tion (except very unlikely Dirac-like size distributions), these ero- 
sive impacts outnumber by far the low-velocity s\ ~ s 2 encounters 
and are the dominant collision outcome. Only the innermost region 
within ~ 0.5 AU has (Av) low enough to allow planetesimal accre- 
tion. 

Interestingly, these results are relatively robust when varying 
the gas disc density. The main difference is that the outer limit of 
the inner "accretion-safe" region varies between 0.35 and 0.75 AU. 
Only for very gas poor systems do we get a qualitatively different 
behaviour, with the safe region extending up to ~ 1.3 AU, where a 
high-(Av) zone due to pure secular perturbations begins. 

These results are further strengthened by the fact that our sim- 
plified assumption for the gas disc probably results in a slight un- 
derestimate of <Av)(s 1jS2 ) values (see discussion in Sec l2.2t . As a 
consequence, our approach can be regarded as conservative regard- 
ing the extent of the impact velocity increase which is reached, the 
"real" system being probably even more accretion-hostile than in 
the presented results. 

This hostility to kilometre- size planetesimal accretion is of 
course a major threat to planet formation as a whole, of which it 
is a crucial stage. However, before reaching any definitive conclu- 
sions, one has to see if there might be some ways for planetary 
formation to get round this hurdle. In the next 3 subsections we 
critically examine three such possibilities. 

4.1 Accretion from large initial planetesimals? 

We have seen that, apart from an unrealistic gas-free case, the only 
way to get an accretion friendly environment at 1 AU is to assume 
large, > 30 km planetesimals which can sustain the high velocity 
regime without being preferentially eroded (Fig|7]l. This solution 
to the high velocity problem raises however important issues. One 
of them is how realistic it is to assume a population of large initial 
planetesimals. As already pointed out, the planetesimal formation 
process is still too poorly understood to reach any definitive conclu- 
sion regarding what an initial planetesimal population looks like, so 
that values in the 30-50 km range cannot be ruled out a priori. Nev- 
ertheless, one has to ask the question of what was before these large 
planetesimals. Should these objects arise from the progressive ac- 
cumulation of grains, rocks and pebble s in a turbulent nebula (as i n 
the standard coagulation model, e.g. lDullemond & Dom inik 2005), 
then there is no clear cut separation between the planetesimal- 
formation phase and the subsequent planetesimal-to-embryo phase 
|fl In this case, the stage where the biggest planetesimals are ~50km 
in size should be preceded by a stage where the biggest objects 
are in the 1-10 km range, which leads us back to the accretion- 
inhibiting nominal case. The only solution would be to bypass these 
intermediate accretion-inhibiting stages by directly forming large 
planetesimals. This could in principle be achieved with the alterna- 
tive gravitational instability scenario (e.g lGoldreich & Wardll 19731 : 



5 In this case, the very concept of an "initial" planetesima l population 
might be questioned (see for instance Wetherill & Inaba 2000) 



lYoudin & Chiand2004) . but t his scenario still h as several problems 
to overcome (see for instance |Xrmitaae 20q3, for a brief review). 
Moreover, it is not clear how gravitational instabilities should pro- 
ceed in the highly perturbed environment of a binary. 

It should also be noted that even if at some point this 
accretion-friendly large-planetesimals stage is reached, the way ac- 
cretion will proceed from there would be very different from the 
standard runaway-accretion scheme. Indeed, as can be clearly seen 
in Fig|7] the (si,Sz) parameter space for which accretion can pro- 
ceed unimpeded is restricted to a narrow 5, ~ s 2 band. For most 
impacts on s > 30 km targets, even if relative velocities are below 
v *(j],j 2 )iow (the "blue" zone), they are nevertheless increased rela- 
tive to the unperturbed case. These increased (Av) would strongly 
reduce the gravitational focusing factor which is the key process 
driving runaway accretion. In this case, the gr owth mode could 
be the so-called "type II" runaway identified bv lKortenkamp et"al] 
(2001), corresponding to a slowed-down runaway growth in a high 
(Av) environment. 



4.2 Accretion once gas has dissipated? 

An interesting issue is what happens after gas dispersal. Indeed, 
it is believed that primordial gas discs dissipate, due to photo- 
evaporation a nd/or viscous accretion onto the central star, in a few 
10 6 years (e.g lAlexander et al.l2006h . a timescale which could even 
be shortened in binaries because of the outer truncation of the 
protoplanetary disc. One might thus wonder if, despite of a long 
accretion-hostile period while gas is present, planetesimal accre- 
tion could start anew once the environment becomes gas free and 
differential phasing effects disappear. However, for this to occur, 
there is a big difficulty to overcome. Indeed, when the environment 
becomes gas free, all planetesimals are still on the differentially- 
phased orbits they were forced onto by gas friction. This means that 
the initial dynamical conditions are such as (Av), si S2) are high. Un- 
der the sole influence of the companion star's secular perturbations 
(the only external force left in the gas free system), these initial or- 
bital differences cannot be erased or damped. The size-independent 
orbital phasing which will be forced by secular effects will come 
in addition to the initial orbital differences, which can in this re- 
spect be regarded as an initial free component. This behaviour can 
be easily checked by gas free test simulations. 

Moreover, even if a hypothetical additional mechanism acts 
to damp the initial differential phasing, it has to act quickly in or- 
der to lead to low impact velocities. Indeed, if given enough time, 
secular perturbations alone will eventually induce high (Av> due to 
orbital crossing of neighbouring orbits. The inner limit of ~ 1.3 AU 
for such orbital crossing given in Sec. |3.2| is indeed the one reached 
at t = 10 4 yrs, but this limit is moving inwards with time. Using 
the empirical analytical expressions derived by TMS06 (see Eq.15 
of that paper), we find that the orbital crossing "front" reaches the 
0.5 AU region after ~ 1.5 x 10 5 yrs. 

As a consequence, it seems unlikely that (Av) (si S2) can be sig- 
nificantly reduced in the later gas free system. Collisional erosion, 
started in the gas-rich initial disc, should thus probably continue 
once the gas has dissipated^. 



6 As this paper was being reviewed, a yet unpublishe d study by Xie and 
Zhou (ApJ, submitted, with a prelimimary summary in Xie & Zhou 2008) 
was brought to our attention. It explores the possibility, for the specific 
case of y Cephei and under possibly extreme assumptions for gas dispersal 
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4.3 Outward planet migration and different initial binary 
configuration 

It could be argued that planets could have formed elsewhere and 
later migrated to the r > 0.75AU regions (the maximum extent of 
the accretion-friendly region obtained in our gas-disc profile explo- 
ration). Early migration of the protoplanets by interaction with the 
primordial gas disc mig ht be ruled out since su ch a migration (be 
it of type I, II or III, e.g P apaloizou et alj|2007» would be inwards. 
In the present case, the planets would thus have to come from fur- 
ther out in the disc, i.e, from regions which are even more hostile 
to accretion. Outward drift of planets could however occur for an- 
other type of migration, the one triggered, after gas-disc dispersal, 
by interactions between a multi-planets system and a massive disc 
of remaining planet esimals (as in the so-called Nice-model, e.g. 
iTsiganis et al.ll2005r) . Netherless, it remains to see how this sce- 
nario, which was studied for mutually interacting giant planets at 
radial distances between 5 and 30 AU, would work in the present 
case, ~ 10 times closer to the star and with probably much less 
massive planets. This important issue will be addressed in a forth- 
coming paper (Scholl et al., in preparation). 

An alternative possibility would be that the binary's orbital 
configuration was not the same in the past as it appears today. 
If the presently observed binary system was part of an unstable 
multiple system (triple or more), planets might have formed when 
the binary had a larger semi- major axis and it was less eccen- 
tric (M arzan & Barbierill2007h . Both these conditions could lead 
to a more favourable environment for planetesimal accretion than 
present. In this scenario planets should have formed prior to the on- 
set of the stellar chaotic phase which inevitably ends with the ejec- 
tion of one star, leaving the remaining binary system with smaller 
separation and higher eccentricity. The discussion of the possible 
consequences of a transient triple system stage for a Centauri A 
and B is complicated by the presence of Proxima Cen, a red dwarf 
located roug hly a fifth of a light-year from the AB binary. A re- 
cent paper by Werthe imer & Lau ghlin ( 2006) argues that the bind- 
ing energy of Proxima Centauri relative to the centre of mass of the 
a Centauri binary is indeed negative, even if the star is close to the 
outer border of the Hill's sphere of the AB system in the galactic 
potential. This might be an indication of a violent past of the sys- 
tem when the AB binary was possibly part of a higher multiplicity 
unstable stellar system. This issue exceeds the scope of the present 
paper and should be explored in future studies. 



4.4 Conclusion and Perspectives 

Because of the hypothetical additional effects listed in Sec 14.31 
the presence of planets beyond 0.75 AU cannot be fully ruled out. 
However, we think that our results on planet formation, with the 
present binary configuration, are relatively robust. Our main result 
is that it is very difficult for s < 30 km planetesimals to have ac- 
creting encounters beyond 0.5-0.75 AU from the primary, which 
makes planet formation very unlikely in these regions. The se re- 
sults are in sharp contrast to those o flMarzari &Scholll j2000). who 
found the terrestrial region aournd a Centauri A to be favourable to 
planet esimal accumulation. The main reason for this discrepancy 
is that iMarzari & Scholll {2000) restricted their pioneering study 



to collisions between equal-size bodies, for which gas friction in- 
deed decreases (Av) and favours accretion, thus implicitly neglect- 
ing the dominant effect of eroding S[ t jf 2 impacts. Our r esults 
are also in contrast to the c onclusions of iBarbieri et alj d2002l) and 
lOuintana et id] d200ll2007l) . who regarded planet accretion as been 
possible in the inner < 2.5 AU region. This is because these studies 
focused on the later embryo-to-planet phase, thus implicitly assum- 
ing that the preceding planetesimal-to-embryo phase was success- 
ful. The present study shows that this is probably not the case. This 
confirms that the planetesimals-to-embryos phase is more affected 
by the binary environment than the last stages of planet formation. 

Due to the complexity of the coupling between secular per- 
turbations and gas drag, and their strong dependence on binary or- 
bital elements, the present results cannot be directly extrapolated 
to other systems. Clearly, such detailed studies of planetesimal ac- 
cumulation in other close binary systems, either those known to 
harbour planets (like y Cephei[j or HD41004) or those being po- 
tential candidates, should be undertaken in the close future. In this 
respect, one object of particular inte rest could be a Cen A's com- 
panion, Alpha Centauri B, for which iGuedes et ail d2008l) very re- 
cently explored the possibilities for observational detection of ter- 
restrial planets. 
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